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TECHNICAL PUBLICATION 


EVOLUTIONARY OPTIMIZATION OF A GEOMETRICALLY REFINED TRUSS 

1. INTRODUCTION 


Structural optimization is a field of research that has received much attention for many years. 
Researchers have used various programming techniques to solve structural problems, for example genetic 
algorithms (GAs), gradient methods, and numerical methods as discussed by Benford and Tinker. 1 A tra- 
ditional problem used in structural optimization is the benchmark problem given in figure 1. This problem 
has been solved using gradient methods for optimization, numeric methods, and GAs. All of the above- 
mentioned methods used cross sections of the individual members as the design parameters. There have 
been many research papers focusing on the cross-sectional area sizing of the individual truss members in 
the benchmark truss problem shown in figure 1. Schmitt and Miura demonstrated a solution to this bench- 
mark problem. 2 ’ 3 They developed an advanced automated design procedure for minimum weight design of 
structures, known as ACCESS. Their solution used finite element analysis (LEA), coupled with nonlinear 
techniques. Schmitt and Larshi created an elemental minimum-weight design tool that utilizes approxima- 
tion steps to replace the design variables. 4 The steps used by these authors are temporary deletion of non- 
critical elements, design variable linking, and Taylor series expansion. Venkayya used an energy criteria 
based on static loading of the structure. 5 He describes this method as extremely efficient when compared 
to nonlinear and linear methods. Using an optimality criteria instead of a search technique has created 
the best design of all the authors and with less computation time. Arora and Haug use a design sensitivity 
analysis to find the minimum cross-sectional members in the benchmark problem. 6 The results achieved 
for the benchmark problem by all of these authors is given in table 1. 

Structural optimization techniques fall into four main categories— topology, shape, size, and mate- 
rial relaxation. Topology optimization is the removal and addition of material within a given design space. 
Shape optimization is the changing of certain dimensions of a fixed topology, for example a picture frame 
and a tire have the same topology but different shapes. Size optimization is a change in the spatial dimen- 
sions of the individual elements in the design region. Linally, material relaxation is the property relaxation 
of certain nonspatial parameters, such as density, modulus of elasticity, and porosity. The optimization 
approach used in this research is a combination of topology and shape. 



Figure 1. Ten-bar benchmark mass minimization truss problem. 


Table 1. Comparison chart for optimal total mass for benchmark problem. 





Optimum Cross Section Areas of Bars, (in 2 ) 


Bar 

Number 

Length 
of Bars 
(in) 

NASTRAN 

Schmitt-Miura 

Schmitt-Farchi 

Venkayya 

Haug-Arora 

Hull etal. 

1 

360 

24.37 

24.43 

24.25 

23.4 

23.27 

Example : 

2 

509 

20.818 

21.06 

20.69 

21.08 

21.2 

46 x 23 

3 

360 

30.62 

30.66 

33.42 

30.41 

30.03 

See fig . 10 

4 

509 

8.4155 

8.58 

8.39 

8.69 

7.47 


5 

360 

0.1 

0.1 

0.1 

0.1 

0.1 


6 

360 

0.22981 

0.1 

0.1 

0.13 

0.1 


7 

360 

0.16575 

0.1 

0.1 

0.1 

0.56 


8 

360 

14.997 

14.59 

14.26 

14.9 

15.29 


9 

509 

0.23011 

0.1 

0.1 

0.19 

0.1 


10 

509 

20.44 

21.06 

20.69 

21.08 

21.2 


Total Mass (lb) 

- 

5,078 

5,074 

5,108 

5,053 

5,010 

5,133 


In addition to the cross-sectional design parameter problem, there has also been much work in the 
structural topology design parameter problems. A topology description of a structural design is a means to 
describe the material connectivity of a structure as well as the number of holes as discussed by Howell in 
a prespecified design domain. 7 Different topology examples are depicted in figure 2. A typical structural 
topology design domain consists of a discretized element mesh and certain properties, such as spatial 
dimensions, Young’s modulus, and Poisson’s ratio. The topology approach to designing structures is typi- 
cally regarded as the first level and most wide-ranging type of optimization. Using this method requires 
the definition of a parameterization scheme for the design space. The topology method for structural 
optimization problems is fathered by Mitchell, and more recently by Kirsch. 8 ’ 9 It has since been applied 
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to compliant mechanism optimization by a variety of researchers including Zhou and Rozvany, Frecker et 
al., Joo and Kota, Poulsen, Pederson et al., Bruns and Tortorelli, and Ananthsuresh and Kota to name just a 
few. 10-16 Structural optimization using this method is essentially a consideration of all possible combina- 
tions of the material connectivity, as defined by the design parameters in the design domain, which meets 
the basic design criteria. This procedure brings a multitude of possible designs into view of the designer 
with the optimum material connectivity for a specified task to be determined by the optimization proce- 
dure. The essence of structural topology optimization is the parameterization of the design space chosen. 


Output 

Disp. 



Figure 2. Topology examples including: (a) design region for a compliant 
mechanism with boundary conditions, (b) input and output forces, 
and (c) various topologies of this problem. 17 


Parameterizing the structural design space is essentially dividing it into discretized elements of 
various shapes and connectivities. The most widely used discretized elements are blocks. 18-20 Other com- 
mon elements are frame 21 and honeycomb, 22 and they will not be discussed further here. The research 
demonstrated in this Technical Publication (TP) builds on a previous publication by the lead author on 
compliant mechanism design. 23 The work presented here follows a traditional structural optimization 
format with few changes: (1) objective function definition, (2) analysis of design set, (3) implementation 
of a search or approximation algorithm, and (4) definition of the design parameters for the problem. The 
additional steps included in the traditional design method are: multilevel geometric refinement steps of 
subdivision, solid mechanical analysis using ANSYS®, and a differential genetic algorithm optimization. 

This TP will describe the alternate structural design tool specifically highlighting the following 
steps in the design process: geometric refinement through subdivision, solid mechanics analysis using 
ANSYS, and design optimization with a differential genetic algorithm. This TP will conclude by present- 
ing an alternate solution to a variation of the traditional benchmark problem for a loading case. 
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2. DESIGN TOOL DESCRIPTION 


The design tool demonstrated in this TP is briefly described here; for a more detailed description 
please refer to Hull and Canfield. 23 The structural optimization approach used in this research is based on 
the continuum model approach commonly used in compliant mechanism design. 

The continuum model approach defines a material region that is a representation of all possible 
solutions to the specific compliant mechanism problem. Certain optimization techniques are borrowed 
from the structural design field to distribute the material in the region into an optimal topology based on 
a predetermined objective variable. The geometry or topology of the truss continua for a material is opti- 
mized based on a given input parameter and objective value. A continuous optimization technique, such as 
gradient-based or genetic programming, is used to find a suitable solution to this problem. A material con- 
tinuum design region with boundary conditions and input and output parameters is presented in figure 3. 
The final optimized solution is a combination of topology (connectivity of the elements), shape (boundary 
of the device or specific members), and size (discretized size of the elements, thickness, width, etc). 



Figure 3. Material continua design domain and problem parameters for a compliant 
mechanism using the continuum model approach. 23 


The design tool presented here follows the general outline as the continuum model approach with 
three exceptions: (1) addition of geometric smoothing, (2) solid mechanics analysis through FEA soft- 
ware, and (3) differential genetic algorithm optimization. 

2.1 Geometric Subdivision 

Multilevel geometric subdivision of the populated design region is applied in this proposed optimi- 
zation method to remove any singularity points and high-stress regions. A singularity point is a geometric 
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region in which two squares are connected by a single point. These anomalies are impossible to manu- 
facture and introduce numerical instabilities in the finite element evaluation of a given design. 24 The high 
stress regions in the geometry of the problem are commonly created by sharp interior comers; rounding 
these regions through subdivision overcomes this difficulty. An example of second-level geometric subdi- 
vision is demonstrated in figure 4. 




Figure 4. Two different examples of second-level subdivisions. It is clear how 
the sharp edges and singular points are reclaimed with material. 


Subdivision is a method used to define a curve or surface through a sequence of successive refine- 
ments. The basic concept is traceable to the early 1940s, however the first concept of generating a smooth 
curve from a control polygon was introduced by Chaikin. 25 Many researchers since that time have pre- 
sented subdivision schemes for modeling smooth surfaces of polygons— Doo and Sabin 26 and Catmull 
and Clark 27 both presented the foremost authoritative schemes. These papers presented methods to extend 
bicubic and biquadratic B-splines to arbitrary meshes. Many different subdivision schemes have fol- 
lowed. 


The geometric refinement step smoothes the design boundary from sharp right-angle definitions 
to a more rounded shape. Imagine sanding an exterior corner of a wooden block and filling the interior 
comer created by two wooden blocks, thus removing all sharp exterior and interior angles. The subdivi- 
sion method used in this research is theoretically described as a series of control points that are spatially 
positioned within the boundaries of the design region. The geometric refinement method of subdivision 
creates a new set of control points that will ultimately be calculated from the initial set through a set of 
averaging equations. These subdivision faces, created from the new control points, are defined as the limit 
of an infinite refinement process. The refinement process begins with the original set of control points 
defined over the design region. The additional control points are calculated at successive iterations using 
the relation 

p./ +/ = Sp 7 . (1) 
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Where p is given as the original set of control points, S is the local subdivision matrix and j, j+1 
are the current and newly subdivided levels. Thus, the subdivision matrix S in a simplistic sense defines 
the subdivision process. The S matrix is simply a weightings matrix that governs the positional influence 
of the surrounding points on the newly calculated point. There are many subdivision schemes available; 
this one used here is a Doo-Sabin Midedge scheme. An initial block topology formulation, followed by a 
second-level subdivision application, is demonstrated in figure 4. 


2.2 Finite Element Analysis 

The objective function is evaluated as the solution to the elastic boundary-value problem using 
commercial FEA software. This work employed ANSYS 8.0 for all FEA. The FEA software is called as 
part of the optimization routine and information is sent on a specific design to be modeled in solid model 
form. Pertinent loading and boundary conditions are also sent. In addition, specification on meshing can 
be passed on as well. Automatic variable meshing using a two-dimensional, eight-node structural solid 
element is specified. This element was chosen due to the accurate results for irregular shapes and X and 
Y translation freedoms. Upon completion of the preprocessing and postprocessing stages, ANSYS returns 
solution information that includes nodal coordinates, primary stresses, and nodal displacements. An elas- 
tostatic analysis is used within the ANSYS program to determine the deflections and strains. Using this 
information, the fitness function value is determined for that particular design. While the analysis of a 
single member of the population is performed very quickly, the fitness evaluation of a population— rather 
than a single point in the design space— is required by the GA and forms a primary source of the compu- 
tational costs of this approach. However, the problem lends itself readily to being processed in a parallel 
fashion, where multiple seats of the FEA software are available. Therefore, the nature of this problem 
could reduce some of the negative factors, such as the long computation time, by providing a simple 
means to easily harness a greater amount of computational resources. The evolutionary search method 
used in this research is presented in the section 2.3. 

2.3 Evolutionary Computation 

An evolutionary computation (EC) search method is used in this design method to obtain a better 
solution to a difficult problem. For example, a design region consisting of 50 x-divisions and 50 y-divi- 
sions of 50 results in a search region of 2( 50 * 50 ) = 2 2 ’ 500 , this infinite search space holds many local mini- 
mums. A randomly guided search method such as a G A predictably finds designs in the global minimum 
region. A presentation of the method used in this research is given in section 2.4. 

Evolutionary computation is the field of study devoted to the design, development, and analysis of 
problem solvers based on simulated evolution. 28 ECs have been used successfully to solve a wide range of 
complex optimization and machine-learning problems for which traditional problem solvers yield unsat- 
isfactory results. ECs differ from most traditional problem solvers in that they operate on a population of 
candidate solutions (CSs) rather than on just one CS. 

Figure 5 provides a pseudocode example of an EC. Initially, a population of individuals in the form 
of CSs is randomly generated and evaluated using a user-specified evaluation function. The evaluation 
function assigns each individual a fitness that is representative of its goodness. After the initial population 
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has been created and evaluated, the EC iteratively refines the population by: (1) selecting parents based on 
their fitness, (2) allowing the selected parents to create offspring through crossover (sexual reproduction) 
and/or mutation (asexual reproduction), (3) evaluating the offspring, and (4) determining which individu- 
als of the current population and the set of offspring get to survive to the next generation. 


Procedure EC() 

{ 

t = 0; 

Initialize Pop(t); 

Evaluate Pop(t); 

While (Not Done) 

{ 

Parents(t) = Select_Parents(Pop(t)); 
Offspring(t) = Procreate(Parents (t)); 
Evaluate(Offspring(t)); 


Figure 5. Evolutionary process of an EC. 


Two well-known EC paradigms that are compared in this TP are: (1) GAs 28 and differential evolu- 
tion (DE). 29 - 32 


2.4 Genetic Algorithms 

Figure 6 provides a pseudocode version of an elitist generational GA. 29 In an elitist generational 
GA, an initial population of IP(t)l individuals is randomly generated and each member of the population 
is evaluated. During each iteration of the evolutionary cycle IP(t)l individuals are selected to be parents 
(with replacement) and these parents are grouped into IP(t)l/2 pairs where each pair creates a total of two 
offspring by using uniform crossover. In uniform crossover, for each gene of parents and offspring, one 
value of a parent is given to one offspring with probability 0.5 while value of the other parent is given to 
the second offspring. Figure 7 provides an example of uniform crossover. After crossover, the offspring 
are mutated slightly (to add a little diversity to the population) and then evaluated. The next generation of 
individuals, for the next iteration, is composed of the best individual discovered by the search so far plus 
the best IP(t)l-l offspring. 
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Procedure GA{ 
t = 0; 

Initialize P(t); 

Evaluate P(t); 

While (Not Done){ 

Parents(t) = Select_Parents(IP(t)l from P(t) with 

Replacement); 

Offspring(t) = Procreate(Parents(t)); 
Evaluate(Offspring(t)); 

P(t+1) = (Bcst lndividual u (Offspring(t)-l}}; 
t = t + 1 ; 


Figure 6. Evolutionary process of an elitist generational GA. 


Parent 1: 

xxxxxxx 

Parent 2: 

Y YY YY Y Y 

Offspring 1: 

X YX YYX Y 

Offspring 2: 

YX YXX YX 


Figure 7. Example of uniform crossover. 


2.5 Differential Evolution 

Differential evolution is another well-known, successful EC paradigm. 29-32 DE is used in this TP 
to improve the overall search method by reducing the search time. The DE used in this TP is similar to the 
elitist GA described earlier with the difference being that instead of selecting two individuals to be parents 
of two offspring, three parents are randomly generated and are used to create one offspring. The offspring 
created are a composition of differences of two of the three parents as follows: Offspring^ = Parentl k 
+ yrnd(0,l)(Parent2 k - Parent3k), where i represents the zth offspring and where the parents are vectors 
and the offspring is simply the resultant vector formed by taking a component, k (for all k), from the first 
parent vector and adding to it the difference of the corresponding components in the second and third vec- 
tors multiplied by a random number in the interval [0..] where gamma is a user-specified parameter. Once 
IP(t)l offspring have been created, mutated slightly, and evaluated, the next generation of individuals— for 
the next iteration— is made up of the best individuals discovered by the search so far, plus the best IP(t)l-l 
offspring. Figure 8 provides a pseudocode version of differential evolutionary algorithm (DEA). 
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Procedure DEA{ 
t = 0; 

Initialize P(t); 

Evaluate P(t); 

While (Not Done){ 

Parents(t) = Select_Parents(3IP(t)l from P(t) with 
Replacement); 

Offspring(t) = Procreate(Parents(t)); // 3 parents per child 
Evaluate(Offspring(t)); 

P(t+1) = { Bcst_Individual u (Offspring(t)-l}}; 
t = t + 1 ; 

} 

} 


Figure 8. Evolutionary process of a DEA. 


9 


3. STRUCTURAL OPTIMIZATION PROBLEM FORMULATION 


3.1 Design Space and Problem Parameterization 

The series of control points, which are used to discretize the design space, are initially equally 
spaced over the design space. This may result in locations similar to the nodes of elements commonly 
used to represent the design space or allowable region for the compliant mechanism (CM). 33 The control 
points contain information including position coordinates, material properties, and boundary classifica- 
tion. Based on this control point parameterization the optimization problem can be described as 

min /( p) . (2) 

subject to the control point classifications as defined by a two-dimensional or a three-dimensional 
parameterization scheme. 


3.2 Minimum Spanning Tree Connectivity 

There are many techniques to determine the connectivity of the blocks. The primary author pre- 
viously used a reachability technique from graph theory in Hull and Canfield. 23 While this method is 
effective and very accurate, the computational time required is significant. An adjacency matrix must be 
calculated for every design and inverted. The current method for determining the connectivity for a given 
design is called a minimum spanning tree. 


The minimum spanning tree formulation is typically applied to computer and phone networks to 
determine the shortest path to connect all of the services in a predetermined network. Also this method is 
used in a more discrete way by solving the traveling salesman problem. The minimum spanning tree is 
used here to simply check connectivity of the blocks in a given design. 


The applied method using the minimum spanning tree is briefly described here. Let X be a vector 
described by the coordinates of H. Then let m be the shortest edge that connects X to H. Therefore m is cat- 
egorized as part of the minimum spanning tree. Subsequently the entire problem is approached by finding 
the shortest distance between two points. When applied to the topology connectivity problem a minimum 
spanning tree is determined to exist or not between necessary blocks, for example, between boundary 
condition and force application blocks. 
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4. BENCHMARK EXAMPLE PROBLEM 


The example problem from the literature will be used to compare the alternate solid model design 
method to the well-known 10-member truss problem. The traditional problem setup that includes geo- 
metric boundaries and force and boundary locations is shown in figure 1. The design region with the 
design parameter blocks to be added and removed is shown in figure 9. Some combination of the design 
parameter blocks will be subdivided and evaluated using FEA. The material is assumed to be aluminum 
(10 + E07 psi), with an allowable stress limit of 25,000 psi on all members. A single-load case is consid- 
ered and is given as Pj = 10,000 lb, P 2 = 10,000 lb. Three example problems are presented in this section 
to demonstrate the capabilities of this tool— a 20x 10, 30x 15, and a 46x23 block discretization. All of 
the following GA structural optimization examples were calculated using a Linux 3.3-GHz, 1-GB RAM 
computer. 


4.1 Example: 200 Block Design 

A total of 2,011x10 design control meshes are used in this example problem, thus a search 
region of 2 (2()x |() ) = 2 20() p OSS jblc solutions. The total computing time was 6 hr. The design region for this 
example problem is given in figure 9. The boundary conditions are applied on the left-hand side and the 
forces applied to the bottom. Here a 150-generation GA search routine is performed. The resulting truss 
design is shown in figure 10. 
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Figure 9. Design region discretization of benchmark truss problem— 5,000 elements. 
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Figure 10. Alternate-deformed optimal topological design for the truss — 200 blocks: 
(a) unsubdivided and (b) subdivided deformed strain plot. 
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4.2 Example: 450 Block Design 


This example design uses 30 x 15 or 450 possible blocks, which gives a design region of 2 450 
(figs. 11-12). A search region of this size took ~12 hr to complete 275 generations using the hardware 
described previously. 



Figure 11. Design region discretization of benchmark truss problem— 450 blocks. 


13 



(b) 



Figure 12. Alternate-optimal topological design for the truss— 450 blocks: 
(a) unsubdivided and (b) subdivided deformed strain plot. 
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4.3 Example: 1,058 Block Design 


A total of 46 x 23 design control meshes are used in this example problem, thus a search region 
of 2( 46x23 ) = 2 1 ’ 058 possible solutions. The total computing time was 47 hr. The design region for this 
example problem is given in figure 13. The boundary conditions are applied on the left-hand side and the 
forces applied to the bottom. Here a 400-generation GA search routine is performed. The resulting truss 
design is shown in figure 14. 



Figure 13. Design region discretization of benchmark truss problem — 1,058 blocks. 


A truss topology design problem was solved in this section. This design is the optimum combina- 
tion of (1) reduction in mass, and (2) minimized deflection and stress for the given design parameters. 
Each of the designs throughout the GA search process had to overcome the singularity conditions associ- 
ated with block formulations, and was accomplished through geometric subdivision. These results show, 
in general, that the truss masses that result from this design tool readily come close to the masses given in 
table 1 by other authors. 

A prototype of the genetically subdivided design truss was machined using a two-axis computer 
numerical control (CNC) machining tool. The material used is lexan for photoelastic properties and the 
design is scaled for viewing purposes, as shown in figure 15. 
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(b) 


Pi 


P 2 


Figure 14. Alternate-optimal topological design for the truss — 1,058 blocks: 
(a) unsubdivided and (b) subdivided deformed strain plot. 
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(a) 




4.4 Example: 5,000 Block Design 

There is no limit on the ability of the code to extend the design parameters to extreme lengths. As 
shown in figure 16, the total allowable blocks are 5,000. This immense parameter set has a search region of 
25,000 t | 1C problem lies in performing the FEA. A design such as this takes in excess of 25 min to perform. 
Also note that an increased solution space will significantly increase the number of design evaluations 
performed through amplified population size and generations. Currently a design of this magnitude is not 
feasible on the computing resources available. Using supercomputing capabilities would still render the 
optimization time in months. 
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Figure 16. Alternate-optimal topological design for the truss— 1,058 blocks: 
(a) unsubdivided and (b) subdivided deformed strain plot. 
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5. DISCUSSION AND CONCLUSIONS 


The preceding examples demonstrate the capability of the evolutionary topology subdivision design 
tool to produce designs on the same order of comparison as those demonstrated in table 1. The topology 
solution to the example problems is a product of the design process presented, specifically engineering 
design, computational elastostatic FEA and an evolutionary search technique. The benchmark-compared 
example 46x23 problem was optimized using =12,000 FEAs within the design region of 2 11058 possible 
designs. The search for the optimal solution to this problem was governed by a randomly guided GA 
search. It is important to realize that the optimality of the presented design is only as good as the number 
of design parameters used. For example, the preceding examples of this tool use fewer design parameters, 
thus the optimality and search regions of these problems is far less. Note the presented topology tool for 
finding optimal truss designs does not include any buckling analysis. 

The design tool presented in this work is ideal for space applications, due the mass savings and 
easy transition to inflatable structures. It is the intent of the author to apply this technology in three dimen- 
sions to address the need for optimal design tools for inflatable structures. As shown in figure 17, the two- 
dimensional design parameterization and subdivision is readily extended to three dimensional. This figure 
is simply a second-level subdivision application to a randomly generated block design domain topology. 



Figure 17. (a) Randomly generated three-dimensional block topology 
and (b) second-level subdivided random topology. 


In conclusion, a genetic design tool is presented using multilevel subdivision and FEA. An exam- 
ple is presented and compared to other examples of a given problem in the literature. While noted that this 
is not a direct comparison, nor is the design better than what is in the literature, but it has been shown that 
this design tool is capable of finding optimal designs. Also, by increasing the number of design parameters 
the objective, mass, and deflection of the design can be minimized further. This tool is developed with the 
intention for applications in space exploration— specifically, thin-film optimal inflatable design. 
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